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Abstract 

We study the quantum fidelity approach to characterize thermal phase transitions. Specifically, 
we focus on the mixed-state fidelity induced by a perturbation in temperature. We consider the 
behavior of fidelity in two types of second-order thermal phase transitions (based on the type of 
non-analiticity of free energy), and we find that usual fidelity criteria for identifying critical points 
is more applicable to the case of A transitions (divergent second derivatives of free energy). Our 
study also reveals limitations of the fidelity approach: sensitivity to high temperature thermal 
fluctuations that wash out information about the transition, and inability of fidelity to distinguish 
between crossovers and proper phase transitions. In spite of these limitations, however, we find 
that fidelity remains a good pre-criterion for testing thermal phase transitions, which we use to 
analyze the non-zero temperature phase diagram of the Lipkin-Meshkov-Glick model. 



PACS numbers: 



I. INTRODUCTION 

Quantum phase transitions (QPTs) the sudden change in the properties of a quantum 
many-body system as a control parameter is varied, have received considerable attention in 
the past decade. Once almost exclusively the domain of condensed matter physics, the 
field of quantum critical phenomena has recently attracted the attention of the quantum 
information community: some quantum entanglement measurements ^ such as concurrence 
[sl, entanglement entropy 4|, and geometric phase [sl can exhibit singular behavior at 
quantum critical points. Thus, they can be used in place of macroscopic thermodynamic 
quantities in classical statistical mechanics - e.g. specific heat and magnetic susceptibility 
- not only to characterize different QPTs, but also to gain insight on the nature of the 
quantum critical behavior. 

Motivated by the sensitivity to perturbations of quantum systems near a critical region, 
one of us and collaborators proposed to use the Loschmidt echo as another quantum 
information probe of QPTs. Based on this work, Zanardi et al further proposed a geometric 
measure: the quantum fidelity (the overlap) between two ground states corresponding to 
slightly different values of the controlling parameters. A flurry of work ensued 9|, showing 
that, despite its simplicity, quantum fidelity does indeed capture the dramatic changes in the 
structure of the ground state at a quantum critical point. In particular, it has been observed 
that for second order QPTs fidelity presents a minimum 

became the standard criterion for detecting quantum criticality with fidelity. Though fidelity 
is used to study QPTs at zero temperature, its finite-temperature (thermal state) extension 
has also been considered [l^. The motivation behind this approach is similar to that for 
QPTs: The proximity to criticality must be reflected in the geometric distance between two 
states separated by a small perturbation (either in temperature or in an external parameter). 



The fidelity of mixed-states Uj at finite temperature also gives useful information about 
the zero-temperature phase diagram Studies of finite temperature transitions using 

this fidelity approach have been reported for specific models, such as the Stoner-Hubbard 
itinerant electron model of magnetism, the BCS model [12I, and also the crossover at finite 
temperature in the one- dimensional transverse Ising model (TIM) [l^ . Nevertheless, we find 
that, in the existing literature, the mechanism for which fidelity can be used to characterize 
thermal phase transitions has not been studied systematically. In this paper we will study the 
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mixed-state fidelity approach in general second-order thermal phase transitions, and explore 
its applicability and limitations. We will focus on non-zero temperature phase diagrams and 
illustrate our arguments with specific examples. Finally, we will also discuss the quantum- 
classical transition of the system when increasing the temperature from a new angle: the 
relation between quantum fidelity and magnetic susceptibility. In the rest of this work, and 
unless explicitly stated, we will use the term fidelity to mean mixed-state fidelity. 

This paper is organized as follows: In Section II, we introduce the finite-temperature 
mixed-state fidelity and study its relation to the analyticity of free energy. In particular, we 
show why fidelity can signal phase transitions by establishing its relationship to specific heat 
and magnetic susceptibility. In Section III we study examples of two types of second-order 
thermal phase transitions - either a divergence or a discontinuity of specific heat at critical 
points-, and discuss the corresponding behavior of mixed-state fidelity. The problems in 
characterizing the second type of transitions with fidelity will be shown. In Section IV 
we discuss other limitations of the fidelity approach: fidelity cannot distinguish between 
phase transitions and crossovers, and at high temperatures thermal fluctuations reduce the 
effectivity of fidelity for picking out critical points. In Sec V, we use the Lipkin-Meshkov- 
Glick model as an example to demonstrate that, despite its limitations, fidelity remains a 
useful pre-criteria for thermal phase transitions due to its simple form. 



II. FINITE-TEMPERATURE FIDELITY AND ITS RELATION TO SPECIFIC 
HEAT AND MAGNETIC SUSCEPTIBILITY 

The mixed-state fidelity of two thermal states with small perturbations in temperature 



and controlling parameter is defined as [10|, 



ill 



T{(3q, Ao; a, Ai) = Try' ^PoPi VPo, (1) 
where the thermal states are written in terms of the Hamiltonian H of the system 

with the partition function 

Z(/3„AJ =Tre-^"^(^'^), (3) 

and where we have perturbations in the Hamiltonian-parameter Ai = Aq + 5\ and in tem- 
perature /So = I/ZcbTq, /?! = 1/A;b(To + 5T). In the following we set Boltzmann's constant 



kB to unity. It can be checked that when both temperatures Tq and Ti decrease to zero, 
the mixed-state fidehty reduces to the ground-state fidehty | {4>gsW\4>gs{^ + ^A)) |, where 
\(pGsW) is the ground state of Hamiltonian -^(A) for a particular value of the controlling 
parameter A . 

When 6\ = 0, we define the temperature fidelity jFg(/?0 5 A, A) = jF(/3o, A; /^i, A), which 
simplifies to 

J^^(/3o, A, A) = , ^ ' (4) 
It can be further proved (see Appendix A) that for small perturbations 5T/T <^ 1 |l4j. Il5l| 



^;3(/3o,/5i,A)^e-^''", (5) 

where C„ = —Td^F/dT"^ is the specific heat at constant field obtained from the free energy 
F of the system. 

When 5T = 0, we can define J^x{P, Aq, Ai) = J-'{P, Aq; P, which can be approximated 
as (see Appendix B) 

J'xiP, Ao, Ai) ^ , 2 ) ^ (g) 

^Z{f3,Xo)Z{P,X,) 

The approximation in Eq. ([6]) is due to the fact that in general H{Xq) and H{Xi) do not 
commute with each other, and is valid only for high temperatures such that f3^6X^ <^ 1. 
7rom Eq. ([6]), and using arguments similar to those used for Eq. ([5]), it can be shown that 
isl (see Appendix C) 

^A(/5,Ao,Ai)^e-— ^ (7) 



where x = —d'^F/dX? is the susceptibility related to an external field of strength A 
(Ao + Ai)/2. 



We see then from Eos. ([5 



order phase transition 



13 



14 



and ([7]) how the fidelity criterion for detecting a second- 



16| plays out for mixed-state fidelity: The minima of JF 
are associated with the singularities of the specific heat and magnetic susceptibility. More 
generally, as we will see below, JF inherits all non-analyticities of the free energy, be them 
divergences or discontinuities in its second derivatives. Therefore, it is reasonable to expect 
that fidelity can be used to study thermal phase transitions, just like traditional criteria 
based on specific heat or susceptibilities. 
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It is interesting to note that from the above Eqs. ([5]) and ([7]) we can also obtain the so 
called perturbation-independent fidelity susceptibilities [l^ 

= (5/3)^^ ~ -^Cy (8) 

--^-fx (9) 

We would like to emphasize that Eq. and Eq. (Q hold approximately only for high 
temperatures, and are a bad approximation for low temperatures and especially for zero 
temperature, where quantum commutation relations are relevant. We will discuss this point 
in Sec IV. Usually the calculation of J^x is much more difficult than that of jFg due to 
the non-commutativity of H{Xq) and H{Xi). In the following we will focus on the fidelity 
for a perturbation in temperature JF^, and its application to second-order thermal phase 
transitions. 



III. FIDELITY IN SECOND-ORDER THERMAL PHASE TRANSITIONS 

Second order thermal phase transitions are characterized by non-analyticities in second 
derivatives of the free energy (e.g. specific heat, susceptibility) with respect to thermody- 
namic variables (temperature and external magnetic fields). According to standard classifi- 
cation , there are two types of non-analyticity that need to be considered: discontinuities 
and divergences (also known as A transition). For ordering purposes we call the associated 
transitions type A and type B, respectively, shown schematically in Fig. [H 

In the following we will discuss the behavior of fidelity near the critical points associated 
to these two types of thermal phase transitions. 



A. Type A: divergence of second-order derivative of free energy 

In this type of transition the specific heat at the critical point is much larger than that 
at other points, and diverges in the thermodynamic limit. From Eq. ([5]) we know that 
the critical point signaled by the maximum Cy will correspond to a minimum of fidelity J^. 
Thus, for this type of systems the decay of fidelity as a function of the parameters (and 
with a fixed perturbation 6T) can be used to characterize accurately the phase boundaries. 
A good example of this situation is the two dimensional (2D) classical Ising model. This 
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FIG. 1: Schematic diagrams of two types of second order phase transition. Type A corresponds to a 
divergence of second derivative of free energy, while Type B corresponds to a discontinuity (jump) 
of second derivative of free energy. In this example we plot the first and second derivative of free 
energy that correspond to entropy and specific heat respectively. Second-order phase transitions 
of Type A are also called A transitions. 

system is described by the Hamiltonian H = —JJ2{ij)^i^j + where means 

sum over nearest neighbor sites, Si = ±1, and A is the external magnetic field. Onsager's 



famous solution 



gives the partition function for zero external magnetic field (A = 0), 



f NT 

Z{p, A = 0) = exp <^ iV In [2 cosh(2/3 J)] + — / d(f)\n 

27r Jq 



(10) 



where K = 2 sinh (2/5J)/[cosh(2/5J)]^. By inserting this into Eq. (jl]), we can obtain the 
fidelity ^/siPo, Pi, ^)\x=o for the 2D classical Ising model (see Fig. 2a). 
The minimum of fidelity agrees well with the analytical result for the critical temperature 
Tc ~ 2. 27 J (see for comparison the specific heat on the right panel). Since fidelity decays 
only on the critical lines, we conclude then that its minimum is a good indicator of criticality 
in this type of transitions. Nevertheless, we would like to point out that the decay of fidelity 
at the critical points becomes less drastic for higher temperatures. This is because thermal 
fluctuations tend to wash out the information about phase transitions encoded in the fidelity, 

6 




0.5 1 1.5 2 0.5 1 1.5 2 

/ J 



FIG. 2: Fidelity (left) and specific heat (right) of the 2D Ising model at zero external field. The 
critical temperature Tc ~ 2.27 J, indicated by the dashed line, is clearly signaled by the minimum 
of fidelity. Because the specific heat diverges at the critical point, the 2D Ising model corresponds 
to a Type A second-order phase transition. 

an effect we will discuss in more detail in Section IV. 



B. Type B: discontinuity of second-order derivatives of free energy 

A common type of transitions is characterized by a discontinuity or jump of second 
order derivatives of the free energy at the critical point. This is the case for instance in 
systems described by a simple Landau-Guinzburg theory [l^]. In such systems fidelity will 
not present in general a minimum at the critical points, but somewhere else in the phase 
diagram. A good example of these type of transitions is the Dicke model, a collection of 

two-level atoms interacting with a single bosonic mode via a dipole interaction with an 

n 

atom- field coupling strength A [19(. The Hamiltonian of the Dicke model can be written as 

H = ujQj^ + uja^a + -^{a) + a){J^ + J^), (11) 



where a and are annihilation and creation operators of the bosonic mode, J^, J+, and 
J_ are angular momentum operators of the total spin of the system, uj and Uq are the 
natural frequencies of the decoupled system, and A is the spin-boson interaction strength. 
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Hamiltonian (fTTj) exhibits both a second-order thermal phase transition [20|] and a quantum 



phase transition 



21 



22| . which has been studied using ground-state fidehty [13|]. Here we 



will study the phase diagram of the Dicke model at finite temperatures using mixed-state 
fidelity, Eq. (jl]). The exact partition function of the Dicke model under the rotating wave 
approximation (RWA) is 22 1 

N 



drre 



2 cosh 



2uj 



'1 + 



(12) 



From this partition function one can obtain that there is a second-order phase transition for 
A > 1 at a critical temperature 2l|, |22[ | 



1 



— tanh ^ ( 



^0 



(13) 



From this partition function, Eq. ( fT2i) . we obtain the fidelity and the specific heat of the 
Dicke model (Fig. [3]). 

There are three aspects to highlight from the fidelity and specific heat shown in Fig. 3. 
First, in the region where there is a thermal phase transition, A > 1, the minimum of fidelity 
does not coincide with the phase boundary line. This is easily attributable to the absence 
of a divergence in the specific heat, which by means of Eq. ([5]) implies that the minimum 
of fidelity need not be correlated to the transition line. Second, fidelity presents minima 
in the region < A < 1, where the system only has a crossover (as seen from the specific 
heat. Fig. 3b). This is again explained by the relation between fidelity and specific heat, 
Eq. ([5]): all maxima of the specific heat (which may not necessarily imply a thermal phase 
transition) will become minima of the fidelity. We will explore more of this point in the 
next section. Third, even though the specific heat has a visible discontinuity at the critical 
points, fidelity changes rather continuously across the phase boundary, specially at high 
temperatures. Here we find a surprise, since from Eq. ([5]) we would expect fidelity to be 
discontinuous at the critical points too. However, as mentioned before, thermal fiuctuations 
affect fidelity strongly, and this discontinuity is washed out for high temperatures. With 
these three observations combined, we see that fidelity is actually not a good indicator of 
criticality for type B phase transitions: it cannot correctly signal the critical points with its 
minima, and is not reliable with discontinuities. Furthermore, as in the < A < 1 region of 
the Dicke model, fidelity might identify simple crossovers as phase transitions. 
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FIG. 3: Fidelity (left) and specific heat (right) of the Dicke model for lo = ujq. On the left plot, the 
minimum of fidelity for fixed A is indicated with a white solid line. The thermal phase transition 
line (dashed line in both plots) in the superradiant phase (A > 1) is indicated by a discontinuity in 
specific heat and, accordingly, also in fidelity. The decay of fidelity in the normal phase (0 < A < 1) 
is due to a crossover instead of a thermal phase transition. The jump instead of a divergence in 
the specific heat indicates that the thermal phase transition in Dicke model belongs to a Type B 
second-order phase transition. 

IV. CROSSOVERS AND THERMAL FLUCTUATIONS 

As discussed in the introduction, geometrical arguments about fidelity in critical systems 
lead us to expect that fidelity will have a minimum at the critical transition points. We just 
saw that this should be extended at least to identify discontinuities in fidelity with type B 
phase transition points (akin to the behavior of ground state fidelity in first order QPTs). 
In this section we explore the following question: is it possible to use fidelity, a quantum 
information tool, to fully characterize a critical system at non-zero temperature, i.e. by 
properly identifying all transition points of the phase diagram? 
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A. Crossovers vs thermal phase transitions 



The free energy of a system is analytic everywhere in the A — T plane except at phase 
transition points. But, there are many "normal" systems without transitions, i.e. where free 
energy is analytic simply everywhere. Nevertheless, this does not exclude the possibility that 
at some points the specific heat can become very large, e.g. at the so called crossover points 
In fact, type A transitions in finite systems look like crossovers that become divergences 
only at the thermodynamic limit. Because of the relation between fidelity and specific heat, 
Eq. ([5]), we expect that fidelity will also have a minimum at the crossover point. This, in 
principle, can be seen as another feature of fidelity, i.e. that fidelity can also be used to 
characterize crossovers [l^l • However, we are interested in the different problem of detecting 
a phase transition using fidelity. 

Let us consider the example of the ID Transverse Ising Model (TIM) with Hamiltonian 
H = —JJ2iLii'^i'^i+i + '^'^i)- The partition function of the system is 23] 



2^exp 



N 



IT 



dk In 



, /2JVl + A2-2AcosA; 
cosh ; — — 



(14) 
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FIG. 4: Fidelity J-'fs{(3o, (ii, X) (left) and specific heat (right) of ID transverse Ising model with 
coupling J in the A — T plane. The white solid line indicates the minimum of fidelity. 



We show a contour plot of fidelity for the ID TIM in Fig. H] (a similar figure can be found 



in Ref . 



13|). In this figure we see a minimum of fidelity following what appears as a phase 
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transition line. We know, however, that this model does not have phase transitions for 
finite temperature (one way to see this is to map the ID TIM into a classical 2D Ising 
model, where the inverse temperature is the effective finite size in the extra dimension of 
the classical system). Thus, fidelity alone may not be able to distinguish simple crossovers 
from proper thermal phase transitions. In order to make this distinction, we must resort to 
traditional statistical mechanics criteria - like the free energy and its derivatives. We show 
in Fig. 4 the specific heat for the ID TIM in the A — T plane, which clearly does not have a 
divergence or a discontinuity. A simple fidelity approach would also have led us to postulate 
a thermal phase transition for the Dicke model for < A < 1 (see Fig. 3), which we know 
does not exist from the exact solution. Even the applicability of fidelity to study crossovers 
is not clear: the "crossover line" found with the minimum of fidelity line is different from 



Tf. ~ |1 — A| obtained in other discussions 
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B. Fidelity and thermal fluctuations 

With fidelity arising from a quantum information approach, it is natural to question 
its behavior for moderate to high temperatures, where quantum effects - such as non- 
commutation of operators - might be obscured. Indeed, thermal fluctuations can wash out all 
information about phase transitions characterized by ground-state fldelity [l^. We already 
saw in the Dicke model and 2D Ising examples of previous sections that fldelity singularities 
become blurred for high temperatures, while the speciflc heat shows a singularity for all 
temperatures. We refer again to flgures 2, 3 and 4 for comparisons between JF and Cy. We 
see that the minimum of fldelity (or its discontinuity) becomes increasingly less prominent for 
larger temperatures, eventually disappearing from the numerical precision. On the contrary, 
speciflc heat is not influenced by thermal fluctuations and is a robust indicator of criticality 
up to very high temperatures. 

Let us give a heuristic analysis of the influence of thermal fluctuations on mixed-state 
fldelity. The perturbation in temperature ST can be expressed as 

^ ~ T T + ST~ T{T + ST) ^ ^ 

Hence, from Eq. (ISl) we can write fldelity as 

Tfi^ e~^^\ (16) 
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From this equation we can see that when temperature increases, the effect on fidehty of 
the singularity of specific heat at critical points will be attenuated. For example, if 
the singularity in C.^ develops slowly with the size of a system, it might be very difficult 

ity with finite size simulations. It is important to highlight 



to detect it reliably using fide 
that the fidelity susceptibility [l5| Xi3 = C'd/(4/5^) = T'^C^/A will not be affected by thermal 
fluctuations at high temperature. Hence, even though fidelity itself may not be a good 
indicator of thermal phase transitions at high temperature, fidelity susceptibilities seem to 
be robust - although this is just because they are proportional to traditional quantities such 
as specific heat and susceptibility. 

V. FIDELITY IN THE LIPKIN-MESHKOV-GLICK MODEL A CASE STUDY 

For all the limitations we have discussed, fidelity decay or jump at the critical points 
is still a necessary condition for a phase transition to exist. Therefore, in the cases where 
fidelity is easier to compute than traditional observables from statistical mechanics - such 
as magnetic susceptibility and specific heat -, it can certainly be used as a pre-criterion 
to explore the phase diagram of a system for potential phase transitions. In order to test 
the predictive power of fidelity for thermal phase transitions, we used it to study the phase 



diagram of the Lipkin-Meshkov-Glick (LMG) [26| model of globally coupled spins with 
an external magnetic field. The Hamiltonian of the LMG model in units of the coupling 
energy is 

i<j i 

where erf, a = x^y^z are the Pauli matrices of the i-th spin, 7 is an anisotropy parameter, 
and A is an applied external field. We approached this problem without previous knowledge 
of its phase diagram, partly to test the usefulness of fidelity, and partly (to be honest) 
out of ignorance. In order to solve this model numerically we used a large spin S = N/2 
representation, 

where Ja = J2iLi'^i f^^ ^ — ^yVy^ is the total angular momentum operator. This is con- 
venient because Hamiltonian (fTS!) does not mix subspaces with different projection of the 
angular momentum, and one just has to diagonalize matrices of size N x N. 
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Indeed, our fidelity studies detected something that appeared to be a thermal phase 
transition in the A — T diagram (the anisotropy parameter 7 turns out to be not very 
important as we will see shortly), see Fig. O We confirmed the existence of a thermal phase 
transition with further numerical calculations of the specific heat and susceptibility, shown 
in Fig. 6, and by a mean field calculation that we present here (we are not aware of such 
a calculation for finite temperature in the literature). Under a mean-field approximation, 




X[J] 



FIG. 5: Temperature fidelity of LMG model for = 800, and 7 = 0.2. Both the thermal phase 
transitions and quantum phase transition at A = 1,T = are indicated by the discontinuity and 
decay of fidelity respectively. Nevertheless, for the thermal phase transition, the discontinuity of 
fidelity deviates slightly from the phase boundary given by the mean-field result (dashed line), 
although this is still within the deviation expected for a finite size system. 
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FIG. 6: Magnetic susceptibility x of the LMG model for N = 800, 7 = 0.2. The phase boundary 
given by specific heat and susceptibility agrees well with that of mean-field result (shown in dashed 
line). Notice however that the discontinuity in x for small A is less pronounced, given that the 
phase boundary lies at almost constant temperature (the reverse happens for the specific heat C„ 
at low temperatures and A ~ 1). 

the LMG Hamiltonian (1171) can be written as 

+7(^r - M,)(aJ - My) + 7M,(ar + a|) (19) 

i 

where = jj: ^Xlili '^j"^ ; '^^ = ^) 1/ is the magnetization along the a direction. We will see 
that this is also the order parameter of the phase transition of LMG model, in analogy with 
the ID quantum XY model. The quadratic terms cancel out, and the Hamiltonian flTIJl) is 
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reduced to 

H = Y.{-M^at - iMya\ - \at) + -Ml + ^M^ (20) 

i=l 

By now, the Hamiltonian is a sum of decoupled single-spin Hamiltonians that can be diag- 
onalized directly. The two eigenenergies are 



£^ = ±^Ml + + A2 + ]^{Ml + 7X)> (21) 



and their corresponding eigenstates are 

1^^^ ^ (M,.-z7M,)|T) + (g+-A)|i) 
J [Ml + 7^2) + - \Y ' 
,^ ^ ^ (M.-z7M,)|T) + (g_-A)|i) ^ ^ 

'(M2 + 72M2) + - \Y ' 

where |t) and |J.) are eigenstates of o"^. The self-consistent equations for the magnetization 
(order parameter) are 

I I (23) 

z z 
where z = e~^^+ + e~^^' is the partition function of the mean field single-spin Hamiltonian. 
Combining Eqs. (12T1) . (!22|) . and (l23l) we obtain the following two self-consistent equations 

tanh [/3v/A^ + M^ + 7W^] 



tanh[/3VA^ + M^ + 7W^^ 

My = , . , — X 7M, 

For 7 7^ 1, the above two equations have nontrivial solutions only when either = or 
My = 0. The two parameters can be further determined by the condition of minimum free 
energy. At absolute zero, the free energy equals the ground state energy £_. It is not difficult 
to find that when 7 < 1, M^ 7^ 0, My = leads to the minimum ground state energy, while 
when 7 > 1, Mj. = 0, My 7^ leads to the minimum energy (and Mj. = My for 7 = 1). For 
example, when 7 < 1, the self-consistent equation is reduced to 



tanh^i ^A'/2 + A2 



(25) 
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In the X — T plane, the phase boundary can be determined by setting the order parameter 
to be zero = 0. Then, the critical temperature as a function of external magnetic field 
A is 



T, = ^ , , O^A^l (26) 



tanh^^ A 

This mean-field result agrees well with the phase boundary obtained by fidelity (Fig. 5) 
and traditional criteria, such as specific heat C„ and magnetic susceptibility x (see Fig. 6), 
because the coordination number of the LMG model is TV — 1 - i.e. it is is big enough to 
ensure the mean-field approximation is reliable. 

Thus, we have detected a phase transition with fidelity, which we then confirmed through 
an analytical result. To the best of our knowledge, this is the first time these events occur 
in this order, and lends support to our discussion above that fidelity is a good pre-criterion 
for testing phase boundaries. 

VI. DISCUSSION AND CONCLUSIONS 

In the above Sec. Ill through Sec. V, we discussed the applicability of JF^ to characterize 
thermal phase transitions, and indicated many of its limitations. Here we would like to 
further consider the applications of J^x and the "quantum" (zero temperature) to "classical" 
(nonzero temperature) transition of the system |27|]. For high temperatures, statistical 
fluctuations dwarf quantum ones, and the importance of uncertainty relations for the ap- 
proximation in Eq. (161) decreases. In this case, the fidelity becomes a function of y and 
the phase transition is classical |27|]. This means that with the increase of temperature, 
the fidelity criteria JF^ for QPT becomes equivalent to the susceptibility criteria for thermal 
phase transitions. Nevertheless, at low temperature, especially at zero temperature, the two 
criteria jFx and x differ dramatically due to the quantum and classical nature of the phase 
transitions. This heuristic analysis agrees well with the result of Refs. I^], that with the 
increase of the temperature the phase transition changes from "quantum" to "classical" . 

In summary, fidelity is a good tool to investigate quantum phase transitions, and has 
been extensively studied. However, when extending to finite-temperature thermal phase 
transition, it faces many limitations: 1) Fidelity decay occurs at both thermal phase tran- 
sitions points and crossover lines, and fidelity alone can not distinguish between them. For 
this, we must fall back on traditional criteria such as the free energy and its derivatives. 
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2) For second-order phase transitions with a divergence in second derivatives of free energy 
(type A), drastic fidehty decay only occurs at critical points, and critical lines can be reliably 
indentified . However, for type B transitions - with a discontinuity instead of a divergence -, 
fidelity decay occurs at many places besides critical points, and maximum decay of fidelity 
may not correspond to phase transition points. Fidelity itself might show a discontinuity, 
but it is easily visible only for low temperatures. Hence, the standard fidelity-criterion for 
second-order thermal phase transitions is more applicable to type A than to type B ther- 
mal phase transitions. 3) In general, the fidelity approach is applicable to low temperature 
thermal phase transitions only. When the critical temperature is very high, fidelity may fail 
to signal the transition because thermal fluctuations wash out all the relevant information 
encoded in fidelity. In comparison, traditional criteria based on free energy are not affected 
by thermal fluctuations and are good for any temperature. 

Before concluding this paper, we would like to point out that, despite its limitations 
for finite-temperature transitions, fidelity can still be a very useful pre- criterion to detect 
thermal phase transitions, especially in systems where we have no prior knowledge about 
its order parameter and symmetries, or even topological thermal phase transitions without 
an order parameter and symmetry breaking. Because of its simple form, we can plot the 
fidelity of the system and then exclude the possibility of thermal phase transitions regimes 
without fidelity singularities. Afterwards, we can focus on suspect areas using free energy 
and traditional criteria to distinguish crossovers from thermal phase transitions. 
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APPENDIX A: SPECIFIC HEAT AND MIXED-STATE FIDELITY 

We will see here the relation between fidelity with a temperature perturbation and specific 
heat. From the standard definition Cy = —Td'^F/dT'^, where F is the free energy, and for 
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a sufficiently small perturbation 5T /T ^ 1, we can approximate 
'F{T + 5T/2) + F{T - 6T/2) - 2F(T)" 



a(T) 



T 



2T2 



In 



(5T/2)2 
Z(T) 



2T, Z(T + 5T/2) 
In ■ 



(Al) 



(5T/2)2 ^Z(T + 5T/2)Z(T-5T/2) ST Z{T-6T/2)' 
Where we have used F = -T\nZ{T). Now, multiplying by 6(3^/(3^ = 6T^/{T + 5Tf, and 
keeping the lowest order terms in ST/T, 



(A2) 



where Pq = 1/{Tq — ST/2), and Pi = 1/{Tq + 6T/2). We thus obtain the relation between 
temperature fidelity JF^ and specific heat 



7(- /3o+/3i ' 



^Z{Po)Z{Pi 



Cv 



(A3) 



APPENDIX B: NON-COMMUTATIVE DENSITY MATRIX 



We look for a simplification of the perturbation in field fidelity, 



J^xiP, Ao, Al) = TrW ^/p^pl^/p^ 



(Bl) 



where pa = exp {—PH{Xa))/Z. Usually, H{Xq) and H{Xi) do not commute with each other 



However, we can use the Trotter- Suzuki formula 



where 

A2(//(Ao),i/(Ai)) 
Thus, we have 



Z(/5,Ao)Z(/3,Ai) 

\[[H{Xo),H{Xi)],H{Xi)]\\ 
Z{p 



to approximate 

< /33A2(i/(Ao), //(Ai))e'3|l^(^°)ll+^ll^("^)ll 



(B2) 



1 

12 



|[[i/(Ao),i/(Ai)],/7(Ao)]|| . (B3) 



^x{P, Ao, Al 



' 2 



^Z{P,Xo)Z{P,Xi) 
Z{P,Xo) 



(B4) 



y/z(/3,Ao + f)Z(/5,Ao-f) 
which is Eq. (E]). The validity condition Eq. flB2p indicates that at high temperature or small 
perturbation (typically P^6X^ <^ 1), the fidelity criteria J-'x for QPT becomes equivalent to 
susceptibility criteria for thermal phase transition and thus the phase transition changes 
from "quantum" to "classical". 
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APPENDIX C: MAGNETIC SUSCEPTIBILITY AND MIXED-STATE FI- 
DELITY 



Similar to Appendix A, we approximate the magnetic susceptibility 

d'F ^ F(Ao + f) + F(Ao- 
aA2 ~ (5A/2) 



Hence we have 



d'F F(Ao + f) + F(Ao-f)-2F(Ao) 



P{6XY {6\y F(Ao + f ) + F(Ao - f ) - 2F(Ao) 



8 8T (5A/2) 



2 



2 In Z(Ao) - In Z(Ao + f ) - In Z{To - f ) 



; In 



2 

ZiXo) 



^Z{Xo + 6X/2)Z{Xo-6X/2) 
From Appendix B, we obtain Eq. ([7]) 



(C2) 



UP, Ao, AO . , - e-^x. (C3) 

Z(Ao + f)Z(Ao-f) 



[1] S. Sachdev, Quantum Phase Transitions (Cambridge University Press, Cambridge, England, 
1999). 

[2] For a review, see L. Amico, R. Fazio, A. Osterloh, and V. Vedrab Rev. Mod. Phys. 80, 517 
(2008). 

[3] T. J. Osborne and M. A. Nielsen, Phys. Rev. A, 66, 032110 (2002); A. Osterloh, L. Amico, 

G. Falci, and R. Fazio, Nature 416, 608 (2002). 
[4] C. Vidal, J.I. Latorre, E. Rico, and A. Kitaev, Phys. Rev. Lett. 90, 227902 (2003); J.I. Latorre, 

E. Rico, and C. Vidal, Quantum Inf. Comput. 4, 48 (2004); Y. Chen, P. Zanardi, Z. D. Wang 
and F. C. Zhang, New J. Phys. 8, 97 (2006);S. J. Gu, S. S. deng, Y. Q. Li, and H. Q. Lin, 
Phys. Rev. Lett. 93, 086402 (2004). 

[5] A. C. M. Carollo and J. K. Pachos, Phys. Rev. Lett. 95, 157203 (2006); S. L. Zhu, Phys. Rev. 

Lett. 96, 077206 (2006); arXiv: 0803.1914. 
[6] H.T. Quan, Z. Song, X.F. Liu, P. Zanardi, CP. Sun, Phys. Rev. Lett. 96, 140604 (2006). 
[7] Z.P. Karkuszewski, C. Jarzynski, and W.H. Zurek, Phys. Rev. Lett. 89, 170405 (2002); 

F. M. Cucchietti, D.A.R. Dalvit, J.P. Paz and W.H. Zurek, ibid. 91, 210403 (2003); R. 

19 



A.Jalabert and H. M. Pastawski, ibid. 86, 246 (2001); F. M. Cucchietti, S. Fernandaz-Vidal, 
J. P. Paz, Phys. Rev. A 75, 032337 (2007); W. H. Zurek, F. M. Cucchietti, and J, P. Paz, 



quant-ph/0312207 Acta Physica Polonica B, 38, 1685 (2007). 



[8] P. Zanardi and N. Paunkovic, Phys. Rev. E 74, 031123 (2006); 

[9] S. Chen, L. Wang, S. -J. Gu and Y. Wang, Phys. Rev. E 76, 061108 (2007); L. Campos Venuti 
and P. Zanardi, Phys. Rev. Lett. 99, 095701 (2007); P. Zanardi, M. Cozzini and P. Giorda, 
J. Stat. Mech. L02002 (2007); M. Cozzini, P. Giorda and P. Zanardi, Phys. Rev. B 75, 014439 
(2007); M. Cozzini, R. lonicioiu and P. Zanardi, Phys. Rev. B 76, 104420 (2007); P. Buonsante 
and A. Vezzani, Phys. Rev. Lett. 98, 110601 (2007); A. Hamma, W. Zhang, S. Haas and D.A. 
Lidar, larXiv:0705.0026l H. -Q. Zhou, J. -H. Zhao and B.Li, larXiv:0704.2940l H.-Q. Zhou, 
larXiv:0704.2945l H. -Q. Zhou, J. -H. Zhao, H. -L. Wang and B. Li, ' arXiv:071 1.46511 M. -F. 
Yang, Phys. Rev. B 76, 180403(R) (2007); Y. -C. Tzeng and M. -F. Yang, Phys. Rev. A 77, 
012311 (2008); S. Chen, L. Wang, Y. Hao and Y. Wang, arXiv: 0801.0020; S. -J. Gu, H. -M. 
Kwok, W. -Q. Ning and H. -Q. Lin, larXiv:0706.2495l 

[10] P. Zanardi, H. T. Quan, X. G. Wang, and C. P. Sun, Phys. Rev. A 75, 032109 (2007). 

[11] A. Uhhnann, Rep. Math. Phys. 9 273 (1976); R. Jozsa, J. Mod. Opt. 41 2315 (1994). 

[12] N. Paunkovic, and V. R. Vieira, Phys. Rev. B 77, 011129 (2008). 

[13] P. Zanardi, L. C. Venuti, and P. Giorda, Phys. Rev. A 76, 062318 (2007). 

[14] P. Zanardi, P. Giorda, and M. Cozzini, Phys. Rev. Lett. 99, 100603 (2007). 

[15] W. L. You, Y. W. Li, and S. J. Gu, Phys. Rev. E 76, 022101 (2007); Min-Fong Yang, Phys. 
Rev. B 76, 180403(R) (2008); Y. C. Tzeng, H. H. Huang, Y. C. Chen, and M. F. Yang, arXiv: 
0804.0537 vl. 

[16] P. Zanardi, H. T. Quan, Xiaoguang Wang, and C. P. Sun, Phys. Rev. A 75, 032109 (2007). 
[17] L. E. Reichl, A Modern Course in Statistical Physics (John Wiley Sons, New York, 1998). 
[18] K. Huang, Statistical mechanics, (John Wiley, New York, 1987.) 
[19] R. H. Dicke, Phys. Rev. 93, 99 (1954). 

[20] C. Emary, and T. Brandes, Phys. Rev. E 67, 066203 (2003). 
[21] K. Hepp, and E. Lieb, Ann. Phys., 76, 360 (1973). 

[22] Y. K. Wang, and F. T. Hioe, Phys. Rev. A 7, 831 (1973); M. Tavis, and F. W. Cummings, 

Phys. Rev. 170, 379 (1968). 
[23] S. Katsura, Phys. Rev. 127, 1508 (1962); P. Pfeuty, Ann. Phys., 57, 79 (1970). 



20 



[24] A. Kopp and S. Chakravarty, Nature Phys. 1, 53 (2005); P. Gegenwart, Q. Si and F. Steglich, 
Nature Phys. 4, 186 (2008); P. Coleman and A.J. Schofield Nature 433, 226 (2005); L. Amico, 
F. Baroni, A. Fubini, D. Patau, V. Tognetti, and Paola Verrucchi, Phys. Rev. A 74, 022322 
(2006); L. Amico and D. Patane, Europhys. Lett. 77, 17001 (2007); A. Cuccoh, A. Taiti, R. 
Vaia, and P. Verrucchi, Phys. Rev. B 76, 064405 (2007). 

[25] M. A. Continentino, quantum scaling in many-body systems (World Scientific, Singapore, 
2001). 

[26] H. J. Lipkin, N. Meshkov, and A. J. Click, Nucl. Phys. 62, 188 (1965); N. Meshkov, A. J. 

Chck, and H. J. Lipkin, Nucl. Phys. 62, 199 (1965); A. J. Chck, H. J. Lipkin, and N. Meshkov, 

Nucl. Phys. 62, 211 (1965). 
[27] M. Vojta, Rep. Prog. Phys., 66, 2069 (2003). 
[28] M. Suzuki, J. Math. Phys. 26, 601 (1985). 



21 



